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Abstract. We have obtained deep spectroscopic observations of several nearby gamma-ray burst (GRB) host 
galaxies revealing for the first time the presence of Wolf-Rayet (WR) stars and numerous O stars located in rich 
and compact clusters or star forming regions. Surprisingly, high spatial resolution imaging shows that the GRBs 
and the associated supernovae did not occur in these regions, but several hundreds of parsec away. Considering 
various scenarios for GRB progenitors, we do not find any simple explanation of why they should be preferentially 
born in regions with low stellar densities. All the examined GRBs and associated SNe have occurred 400 to 800 pc 
from very high density stellar environments including large numbers of WR stars. Such distances can be travelled 
through at velocities of 100 km s~^ or larger, assuming the travel time to be the typical life time of WR stars. 
It leads us to suggest that GRB progenitors may be runaway massive stars ejected from compact massive star 
clusters. The ejection from such super star clusters may lead to a spin-up of these stars, producing the loss of the 
hydrogen and/or helium envelopes leading to the origin of the type Ibc supernovae associated with GRBs. If this 
scenario applies toed text/Sc all GRBs, it provides a natural explanation of the very small fraction of massive 
stars that emit a GRB at the end of their life. An alternative to this scenario could be a binary origin for GRBs, 
but this still requires an explanation of why it would preferentially occur in low stellar density regions. 

Key words, cosmology: observations, galaxies; individual (GRB980425, GRB020903, GRB031203), galaxies: stellar 
content, galaxies: abundances, Stars: Wolf-rayet 

1. Introduction sociation between GRBs and supernovae (SNe) indicates 

in many cases that the parent SN population of GRBs is 
Gamma-Ray Bursts (GRBs) are believed to trace the formed by peculiar type Ibc SNe which have been pro- 
death of massive, short lived stars, providing the most en- ^^^^^ ^ be SN1988bw-like SNe by Folcv ct al. (20^. 
ergetic events in the Universe. This is further supported by ^hey are characterised by high luminosity peaks {Mb 
the discovery that several long-duration GRBs are, indeed, f^.^^ .^g 5 ^ -17) and a high expansion velocity of the 
associated with the collapse of massive stars to a black ^jg^^-g^ (30 000 km s'^). Nevertheless, other types of SNe 
hoK referred to as the collapsar model JGalama eLalJ (^ore regular SNIbc or even SNIIn) cannot be excluded. 



19981 I Stanek et al. | 1200^ i Hjorth et alj |2003J. In this ^he identity of the progenitors of GRBs has thus still to 

model, a rapidly rotating star undergoing core-collapse ^e addressed. However, while the peculiar SNIbc may be 

produces a jetted GRB along the rotation axis, and blows f^^her confirmed as the parent SN population of GRBs, 

up the entire star in an energetic supernova explosion ^^^ ^^^^ favoured progenitors for the collapsar model are 

IJMacFadven &: Wooslevlll999HKlose et alJl2004l) . The as- ^-jj^ Wolf-Rayet (WR) stars 

Send offprint requests to: francois.hammer@obspm.fr JHirschi et al] (|2005J) have investigated WR stars and 
* Based on VLT/FORS2 observations collected at the have shown that stellar evolution models including de- 
European Southern Observatory, Paranal, Chile, programme tailed effects of rotation predict the conditions for GRB 
No. 073.B-0482(A). production via coUapsars, namely the conditions for black 
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hole formation, the loss of the hydrogen-rich envelope 
(such as in SNIbc) and enough angular momentum to form 
an accretion disk around the black hole. Furthermore, if 
only stars of the particular WO subtype - thought to be 
intimately related to SNIc - are considered, the GRB pro- 
duction rate can fairly well be reproduced. In their model, 
GRBs are predicted to occur only over a limited metal- 
licity interval at subsolar values (typically at Z^mc to 
ZhMC, i-e. Z ~ 0.2 — 0.4 Zq). However, when magnetic 
fields are taken into account in thes e models, it may we ll 
be more difficult to produce GRBs llPetrovic et alJl2005|) . 
To circumvent t h is pro blem IWooslev fc Heeen l)2006|) and 
lYoon fc LangeJ l)2005|) suggest rapid rotators at low 
metallicity (typically Z <, 0.05 Zq) as GRB progeni- 
tors, since in these cases a nearly chemically homoge- 
neous evolution and a low stellar mass loss can produce 
the right conditions (high specific angular momentum 
stars with no hydrogen envelope) for the coUapsar model. 
Establishing a precise metallicity limit, below which this 
scenario may work, is however difficult due to uncertain- 
ties i n stellar mass loss rat es and initial stellar rotation 
rates l Yoon fc Langeill2005|) . To estabhsh observationally 
the high mass of long-duration GRB progenitors and the 
collapsar model, and to constrain scenarios like those just 
mentioned, it is imperative to confirm the presence of WR 
stars in the region of the GRB, which has never been 
achieved until now, and to better determine progenitor 
metallicities. 

The first com plete optical study of z < 1 host galaxies 
has been done bv lLe Floc'h et alJ (2003). It revealed that 
GRBs occur in galaxies with low luminosities and blu e 
colors (see also IPruchter et al.lll999t ISokolov et alJl200ll) . 
In addition, there is growing evidence that the majority 
of GRB host galaxies are Lya emitters wit h star forma- 
tion rates (SFRs) between 1 to II Mo/yr l)Fvnbo et alJ 
l2003H.Takobsson et alJl2005l) . All this is indicative of low 
metallicity environments, as confirmed by direct abun- 
dance measurements (Prochaska et al. 2004; Sollerman 
et al. 2005; Hammer et al., in preparation). As a result, 
it leads to some controversy, since GRB hosts were for- 
merly believed to be associated with galaxies with strong 
star formation rates a veraging 100 M^T^/y r from radio and 
sub-mm observations (JBerger et al.ll2003ll . It is unclear if 
the radio emission is simply related to star formation or if 
the sub-mm detection could be affected by the lack of spa- 
tial resolution: this may be questioned b y the absence of 
detec tion of most GRB hosts by Spitzer ( Le Floc'h et alJ 
\200(i\ . 

The analysis of the overall spectral energy distribu- 
tion (SED) shows that the GRB hosts have high specific 
star formation rates (i.e. high SFRs with respect to their 
luminosity) and younger stellar populations than an en - 
semble of all field galaxies (e.g. iGhristensen et al J 120041) . 
However, as these authors note, the derived ages of 50—200 
Myr seem to indicate that GRB hosts are not significantly 
younger than starburst galaxies at similar redshifts. These 
"old" ages measured could be the result of large aperture 
effects or composite stellar populations, i.e. the young pop- 



ulation from which the GRB descends might be diluted by 
older stars. 

Further understanding of the GRB production mech- 
anism requires us to study their environments in de- 
tail. In this paper we present deep spectroscopic observa- 
tions of GRB hosts, and focus our analysis on the four 
most nearby host galaxies: GRB980425 at z = 0.008, 
GRB03I203 at z = 0.1055, GRB030329 at z = 0.169 and 
GRB020903 at z = 0.25. They are bright enough to test 
for the presence of WR stars, as well as to robustly estab- 
lish metal abundances through a direct electron tempera- 
ture measurement. Moreover, they are found to be associ- 
ated with a Type Ibc SN (GRB980425/SNI998bw: Patat 
et al. 2001; GRB03I203/SN20031w: Malesani et al. 2004; 
GRB030329/SN2003dh: Hjorth et al. 2003; GRB020903: 
Soderberg et al. 2005), and when using HST imaging, the 
SN location can be identified with some accuracy within 
the host galaxy. 

Our VLT observations are described in Sect. 2. 
The properties of the environments of GRB980425 and 
GRB020903 as deduced from spectroscopy and imaging 
are discussed in Sect. 3. Based on this data we suggest in 
Sect. 4 a new scenario for GRB progenitor stars. In Sect. 5 
we summarise our result and discuss possible implications. 
Throughout the paper we adopt the ACDM cosmological 
model {Ho=70 km s'^ Mpc"\ nM^O.3 and f^A = 0.7). 

2. Observations and measurements 

Spectroscopic VLT/F0RS2 observations were done in vis- 
itor mode in July 2004 (programme No. 073.B-0482(A)). 
Among our eight targets, two nearby galaxies were ob- 
served, the GRB980425 at z = 0.0085 and the GRB020903 
at z = 0.25 using two differents F0RS2 set-ups (600B 
and 600RI grisms with a resolution R ^ 1300). For the 
GRB980425 host, given the size of the galaxy, the slit was 
placed as it is shown in Fig. I. The GRB980425 host was 
observed with a total exposure time of 1800s and 1500s 
with the 600B and 600RI grisms, respectively, and the 
GRB020903 host 7200s with each grism. 

Data reduction and extraction of optical spectra were 
performed using a set of IRAF procedures developed 
by our team, which allowed us to reconstruct simulta- 
neously the spectra and the sky counts of the objects. 
Spcctrophotometric calibration of each grism was done us- 
ing the same star. Broadband filter images were used to 
compute aperture corrections and check the spectropho- 
tometric calibration. 

Flux measurements were performed using the SPLOT 
package of IRAF. Measurements were performed by two 
of us (F.H. and H.F.) and compared. In the case of 
GRB980245 we also compared our results with those 
performed on spectra for which the background light of 
the ESO I84-G82 galaxy had been removed. All results 
are found to be very similar, which is supported by the 
remarkable consistency of the derived temperatures us- 
ing different chemical species. The electron density was 
derived from the S ii line ratio, the electron temper- 



F. Hammer et al.: GRB 



Table 1. Flux measurements of the GRB host galaxies 







GRB980425 




GRB020903 




SN reg. 


WR reg. 


reg. 4 




on A3727 


4.020 


44.00 


2.680 


6.89 


Nelll A3869 


0.550 


11.30 


0.055 


2.05 


O III A4363 


<0.030 


1.55 


<0.030 


0.54 


NIII A4640 


- 


0.30 


- 


0.0078 


Hell A4686 


<0.052 


0.53 


- 


0.026 


ArlV A4711 


- 


0.34 


- 


0.12 


ArlV A4740 


- 


0.16 


- 


0.20 


H/3 


0.910 


35.40 


0.230 


4.40 


HI A4959 


0.900 


49.70 


0.570 


7.40 


Oni A5007 


2.040 


201.20 


1.100 


33.50 


Nil A5755 


- 


0.17 


- 


- 


1 A6300 


0.400 


1.70 


0.200 


0.65 


SII A6312 


0.040 


0.90 


- 


- 


Nil A6548 


0.380 


4.17 


0.340 


0.25 


Ha 


3.860 


183.90 


1.890 


16.80 


Nil A6583 


0.880 


11.30 


0.520 


0.72 


Sn A6716 


0.600 


10.90 


0.540 


2.39 


Sn A6731 


0.840 


8.60 


0.590 


1.71 


Arlll A7136 


0.190 


7.90 


0.096 


- 


II A7325 


0.195 


3.20 


0.126 


- 


Aril A7751 


0.03 


2.08 


- 


- 



Notes: Measured fluxes before aperture correction in 10 



-16 , 



erg cm A . Typical errors (Icr) are 0.03 for the faint 
lines as estimated from the SPLOT routine. 



ature from O ii 7(3726 + 3729)//(7320 + 7330), O ill 
7(4959 + 5007)//(4363), N ii 7(6548 + 6583)//(5755), and 
S II 7(6716 + 6731)/7(4068 + 4076), when available. These 
estimates were made using the NEBULAR package of 
IRAF, after correcting the lines for extinction. The pack- 
ag e is based on a 5-level atom program and is described 
bv lShaw fc Dufoun l|1994|l . Extinction correction was esti- 
mated using the Ha/H/3 ratio, and verified using higher or- 
der Balmer lines. All have been found in good agreement, 
but one, for which we have corrected flu xes from the un- 
derlyi ng absorption using the method of llzotov fc Thuanl 
l)l999|) based on Ha/H/3 and H7/H/3 ratios. 

3. Properties of the environments of GRB980425 
and GRB020g03 

We now discuss our F0RS2 spectra of the SN1998bw 
remnant and several nearby regions of the GRB980425 
host as well as HST imaging of this galaxy, and 
F0RS2 spectroscopy and HST imaging of the more dis- 
tant GRB020903. A brief description of the hosts of 
GRB030329 and GRB031203 (based on archival data) is 
given for comparison in Sect. 4. 



3.1. GRB980425 host galaxy 

Figure 1 displays the slit position with labels for the re- 
gions around SN1998bw. HST/STIS images (ESO/HST 
archive programmes 8640, 8243, and 8648) of two regions 



are shown, rev ealing 7 sub-areas i n the immediate vicinity 
of SN1998bw llFvnbo et al.ll200nl) . and the "WR region" 
includes a very luminous main component, hardly resolved 
(0.03 arcsec FWHM), corresponding to 5 pc FWHM. 



3.1.1. A bright and compact HI I region with WR stars 

This region lies approximately 0.8 kpc NW from the 
SN1998bw, and its spectrum presents numerous emission 
lines, with a broad range of chemical species (see Fig. 2 
and Table 1). Its density, temperature and oxygen abun- 
dance (see Table 2) are characteristic of a vigorous star 
forming region with Z ~ 0.5 Zq (assuming 12-|-log(0/H) 
= 8.69 for the solar abundance, see AUende Prieto et al., 
2001). Ha and H/3 show very large equivalent widths, 1400 
and 200 A, respectively. After correcting for the aperture 
effect (factor 4.6) and for extinction, we find a luminosity 
of 1.24 X lO^*' erg s^^ in the H/3 l ine, which corresponds t o 
- 2300 equivalent 07V stars f cf. ISchaerer fc VaccJl998|) . 

The spectrum shows a prominent blue WR bump 
around the He Ii4686 line, with many faint emission lines 
(Ariv4711,4740, Niii4640, and Feiii4658, see Fig. 2). The 
absence of carbon hues (Ciii5696 and Civ5808) suggests 
that most of the WR stars are of the WN type. Using the 
luminosity of the blue bump from 4650 to 4686 A (see 
ISchaerer &: Vaccalll998|) . we estimate the number of WN 
stars to be 86 ±20 (Icr error estimated from the contin- 
uum) in that region (assuming WN7 star luminosities). 
It agrees with the measurement of the faint Heii4686 line 
(see Table 1), which corresponds to the emission of 46 to 
92 WN stars. The large Ha and H/3 equivalent widths, the 
presence of WR stars, and fits to the spectral energy dis- 
tribution are all indicative of a very young (~ 1 — 6 Myr, 
depending on the age indicator used) massive star forming 
region. These properti es, as well as the size of the region 
ijO'Connell et al.l[l994|) . show that we are dealing with a 
(possibly multiple?) massive, young, super star cluster. 

We also notice the possible presence of Mgl]4571 in 
emission, a feature c ommonly seen in the nebular phase 
of type Ibc SNe (see IPatat et allboOll: iFolev et al.ll2003(l 
and broad wings in the bottom of the Ha emission line 
revealing gas motion at velocities up to 3000 km s^^. If 
real, these features, combined with the large number of 
WR stars, might indicate the presence of recent super- 
novae in this region. 

While the locati on of SN1998bw is beyond doubt (see 
iGalama et ai]ll999l) . its association with GRB980425 has 
been much debated. The major argument in favour of the 
association of the two events is their remarkable temporal 
coincidence. The fact that other SN of the same type as 
SN1998 bw have been foun d a few weeks after GRB events 
(see e.sf. lFolev et al.l2003|) is also a good argument to sup- 
port the association. Nevertheless, the peculiar properties 
of the WR region near the SN has prompted us to in- 
vestigate whether some supernovae might have occurred 
there few years before our observations. We h ave used 
the HST/STIS imagery fsee iFvnbo et alJl200n|) taken at 
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Table 2. Physical properties of the GRB host galaxies 

GRB980425 GRB020903 

region SN region WR region 4 



N, [S ii]" 


27 


158 


498 


27 


Te [0 II]" 


11939 


11129 


9160 


- 


Te [0 III]" 


< 13966 


11900 


< 18274 


14852 


Te [N II]* 


- 


11249 


- 


- 


12+log(0/H)]^ 


8.25 


8.39 


9.00 


7.97 


N/0 


0.21 


0.043 


- 


0.013 


Av'' 


1.02 


1.51 


0.98 


0.76 



Notes: (a) Electronic density estimated from the S ll flux ratio J(6716)//(6731). (b) Temperature in K estimated from the 
following flux ratios: O II 1(3726 + 3729)//(7320 + 7330), O ill 1(4959 + 5007)// (4363), N ll 1(6548 + 6583)//(5755), 
and S II 1(6716 + 6731)//(4068 + 4076). (c) Metallicity, 12+log(0/H), estimated using the effective temperature method, 
(d) Extinction coefficient, Av (in magnitudes), has been computed using the standard Balmer ratio of H/3 and Ha. 

Fig. 1. SEE ATTACHED JPG FIGURE - VLT image of the ESO-184-G82 galaxy showing the slit location and the 
regions surrounding the SN1998bw remnant. On the right, enlargement of the two area (WR region on top, SN region 
on bottom) based on HST/STIS imagery. The location of the SN is indicated by a circle. 
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Fig. 2. Spectra of the WR region of the GRB980425 host (top, left) and of the GRB020903 host (bottom, left). On 
the right is shown an enlargement around Heii4686 of the spectrum of the WR region of the GRB980425 host. 



three different epochs, namely June and November 2000 
and August 2001. Special care has been taken to derive a 
common astrometric solution for each image and we de- 
rive a 0.03 STIS pixel (or 0.75 mas) accuracy using three 
stars around the WR and SN regions. While the variabil- 
ity of SN1998bw is well detected, we found no indication 
of variability in the WR region. We notice, however, that 
the detection of a faint signal in a very bright area is pre- 
vented by the Poisson noise of the total signal, and so 
we would not be able to see a faint signal such as the 
SN1998bw remnant two years after its maximum, if it was 
embedded in the bright WR region. However, we believe 
that a strong variability similar to that of SN1998bw at 
its peak would have been detected in the WR region, since 
this field has been thoroughly surveyed since May 1998. 



Recall that the SN1998bw luminosity near its maximum 
surpasses the luminosity of the WR region (see Fig. 2 of 
iGalama e t al. 1999). Hence, we conclude that the associ- 
ation of the GRB98G425 with SN1998bw is very probable, 
and that the WR region is a strong star- forming region in- 
cluding a very bright (possibly multiple) super star cluster 
separated by 800 pc from the SN. 

The properties of this cluster are indeed interesting, 
and sufficiently strong to explain reasonably well the de - 
tection at 24 //m reported by iLe Floc'h et alJ l)2006J) . 
Indeed, assuming standard SFR conversion factors, the 
star formation rate derived from the Ha line is 0.34 Mq/jt 
(aperture and extinction corrected), a value in excellent 
agreement with what can be inferred from its IR luminos- 
ity {LiR = 2 X 10^ Lq, see lLe Floc'h et ah 2006) . i.e. SFR 
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= 0.35 Mq/yi. iFiom. its absolute scales the properties of 
the WR region are also similar to those of the embedded 
super star cluster ( "supernebula" ) in the nearby starburst 
NGC 5253. The bolometric luminosity of the dominant 
embedded cluster in NGC 5253 is Lhni ~ (1 - 3) x 10^ Lq 
iJBeck et al.lll996tlVanzi fc Sauvagell2004J) . and it contains 
several thousand O st ars within a small ( pc or even sub- 
pc scale) region (e.g. iTurner et al.ll200J^ . However, the 
extinction in the WR region of the GRB980425 host is 
much smaller. Further comparisons are beyond the scope 
of this paper. 



3.1.2. The region surrounding the SN1998bw remnant 

This region is almost 10 times fainter than the WR region, 
and provides a spectrum with a lower S/N. Nevertheless, 
it shows enough emission lines to derive a full diagnosis 
of its interstellar medium. It shows a more moderate ex- 
tinction and with Z ^ 0.36 Zq a somewhat lower oxy- 
gen abundance than the WR region and region 4 (see 
Table 2). Interestingly, and in contrast to the WR region, 
this region also has a high nitrogen over oxygen abun- 
dance of N/0 = 0.24, almost twice the solar value. While 
such a high N/0 ratio is approximat ely 4 times larger 
than that found in HII regions (cf. e.g. lLiang et alJl2005i 
llzotov et al.ll2006|) . suc h strong N line iri tensities can be 
found in SN remnants IjSmith et al.lll993|) . This N excess 
could be related to the enrichment from the progenitor 
star of SN1998bw, as expected from rotating stars. For 
example, similar N/0 ratios are predicted by the yields of 
~ 10-20 Mp) sta r s in th e rotating stellar evolution models 
of iHirschi et al.l l|2005|) albeit at somewhat larger metal- 
licity than the one observed for this SN region. However, 
dilution with pre-existing ISM will reduce the resulting 
N/0 ratio. 

The SN region shows moderate Ha and H/3 equivalent 
widths of 88 and 16 A, respectively, indicative of ages of 
~ 6-8 Myr adopting our evolutionary synthesis models. 
The aperture (factor 3.27) and extinction corrected H/3 
flux corresponds to only ^ 10 07V equivalent stars. We 
do not detect the presence of WR stars, which cannot only 
be attributed to the relatively low S/N of the spectrum. 
Indeed, Fig. 3 compares the spectrum of the SN region 
with that of the WR region if it was observed at the same 
S/N. The derived upper limit on the luminosity of the 
Heii4686 line (see Table 1) indicates fewer than or 1 WR 
star in that region. This is not surprising given the small 
number of massive stars and the relatively low WR/0 star 
ratio, WR/0 ^ 0.05, observed in the WR region. Recall 
also that HST/STIS (see Fig. 1) is able to resolve the SN 
region in 7 small sub- areas, and that the brightness at 
the precise location of the SN is small comp ared to other 
sub-ar eas. According to the photometry of iFvnbo et alJ 
l)2000|) . in April 2000, the SN contributed to only 13% of 
the V luminosity of the region. Our spectroscopic observa- 
tions were made 3 years after the STIS observations and 
the SN remnant probably has faded away since April 2000: 



it is probable that no massive stars are present in the SN 
sub-region. 



3.2. GRB020903 host galaxy 

The GRB occurred in the outskirts of a small, very ir- 
regular galaxy aX z — 0.25 with Mb{AB) = —19.3 
(see Fig. 4). The slit was c entered on the positio n 
of the reported afterglow (see ISoderberg et alJ l2004al) . 
GRB02 0903 presents ma r iy sirn ilarities with GRB980425. 
Indeed, ISoderberg et alJ 1)200511 have convincingly shown 
that the event was followed by a supernova of the same 
type as SN1998bw. The residual HST image in Fig. 4 (top- 
right) shows the precise location of the SN, which has oc- 
curred at a 0.115 arcsec offset from a compact, unresolved 
region, in the outskirts of the host galaxy. At the distance 
of the host galaxy (z = 0.25), this offset corresponds to 
460 ±100 pc, a value close to the distance between the 
WR and SN regions in the GRB980425 galaxy 

The spectrum at the position of GRB020903 (see 
Fig. 2) reveals an active star forming region, with a strong 
oxygen deficiency oi Z = 0.19 Zq. Ha and H/3 show large 
(rest-frame) equivalent widths, 245 and 40 A, respectively. 
These values are intermediate between those of the SN and 
WR regions of the GRB980425 host. This is not surprising, 
knowing that the slit has likely included several regions 
around the GRB location. After correcting for extinction, 
we find a luminosity of 6.9 x 10^^ erg s~^ in the H/3 line, 
which can be interpreted as produced by ~ 1300 O stars 
(assuming a 07V stellar type) present in the aperture. A 
larger value might be inferred if the aperture correction 
(factor 6) was applied. 

Figure 2 shows a strong signature of a blue WR bump 
and using the Heii4686 emission line, we infer a WNL/0 
ratio ranging from 0.14 to 0.2 (approximately ^ 200 WR 
stars) . The blue bump does not show the same features as 
that of the GRB980425 WR region. This is not so surpris- 
ing after examining the large variety of blue bump features 
from a survey of WR galaxies (Guscva et al. 2000). The 
WR bumps are made of blends of a large variety of emis- 
sion lines and can thus be subject to variations from one 
galaxy to another. The relatively low S/N of our spectrum 
can also alter the appearance of the blue bump. It is more 
difficult to assess the age of the stars near the location 
of the GRB, because the slit likely includes components 
from other regions of the galaxy. Assuming the extinction 
(Ay = 0.76) derived from the ionised gas, the fit of the 
spectral energy distribution reveals the need for an older 
stellar population. 

In summary, the properties of the GRB020903 host 
galaxy seem to be very similar to those of the GRB980425 
host, except that it shows a higher electron temperature 
and a smaller oxygen abundance. As for GRB980425, the 
supernova associated with GRB020903 seems to have oc- 
curred at several hundred parsecs from a bright, relatively 
compact region, which is most likely responsible for the 
numerous WR and O stars seen in our spectrum. 
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Fig. 3. Spectrum of the SN1998bw remnant region (bottom) compared to the scaled spectrum (top) of the WR region 
at the same S/N of the GRB980425 host. It shows that the absence of WR features in the SN1998bw remnant cannot 
be related simply to the lower S/N. 

Fig. 4. SEE ATTACHED JPG FIGURE - Top panels: HST/ACS/F606W image take n 91 (left) and 300 (m iddle) 
days after the GRBG2G903 and showing the associated SN (same type as SN1998bw, see ISoderberg et aObOOSJ) . The 
image on the right shows the residual derived from our extraction, and confirms the result of Soderberg et al. The SN 
(peak identified with the light cross) exploded 460 pc off a luminous region, which is not resolved even at the resolution 
of the ACS. Middle panels: HST/ACS/F606W image taken 8 (left) and 25 (middle) months after the GRB030329 
event and showing the associated SN (right, light cross). Bottom panels: GRB980425 image at z = 0.0086 (left) 
plotted as it would be seen if it was at the redshift of GRB030329 (middle) and of GRB020903 (right). This illustrates 
that the offset between the SN and WR regions in GRB980425 is similar to that observed in the two more distant 
GRBs. 



4. Discussion 



The proximity of the two GRB host galaxies (980425 and 
020903) has allowed us to detail at an unprecedented level 
their spatial and spectral properties. However, these two 
GRB events are amon g the least energetic GRBs (see 
ISoderberg et alJl2004b(l . and in the following discussion 
one should be cautious before extrapolating their proper- 
ties to the numerous cosmological GRBs observed at large 
distances. However, it has been argued that sub-energetic 
bursts are simply events v iewed away from the jet axis. 
Indeed. lGuetta et al.l f2004j have shown that this interpre- 
tation matches the statistics of both low energetic GRBs 
at low red shifts and high ener getic GRBs at cosmological 
distances. iFvnbo et a l. f2004j also argue for such a case 
for the X-ray flash XRF 030723. 



4.1. Inferences on the connection between GRB and 
star formation 

GRBs occur in galaxies that are gene rally sub-L* and 

I 1 1 JL 

even dwarfs (see iLe Floc'h et al.M2003l) . The spectra of 
their hosts include a large variety of emission lines, in- 
cluding most of the Balmer series and Helium lines. This 
is characteristic of a very young burst, which is further 
supported by the fact that few Myr old stars (including 
WR stars) dominate the continuum. Within these small 
galaxies, we find that GRBs (and their associated super- 
novae) do not occur in regions containing a large number 
of massive stars. This is in strong contrast to simple ex- 
pectations from the massive star collapsar model, where 
the GRB progenitor star is thought to represent a rare 
case among a large number of massive stars. 

In both cases studied here, we find that the GRB oc- 
curred at distances of ^ 400 to 800 pc from a compact 
and luminous region, which is just resolved (5 pc FWHM) 
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Fig. 5. Spectrum of the GRB031203 host (2 x 1800 s) taken in September 2004 using the 300V grism with VLT/F0RS2 
(from ESO archive, programme No. 073.D-0255(A)). The red slope of the ene rgy distribution can b e attributed to 
the low Galactic declination {b^^ = 4.7) leading to an extinction Av^ = 3.62 (see lProchaska et al.ll2004J) . A blue bump 
around the Heii4686 line is detected as well as a bump around the CivSSOS line, confirming the presence of WR stars 
in the host galaxy. 



in the case of the WR region in the GRB980425 host. In 
other words, GRBs neither occur directly in massive star 
formi ng regions nor in massive LIRGs (see Lc Floc'h et al, 
120061) . but in regions that show very few or no massive 
stars (e.g. the SN region in the GRB980425 host). This, to- 
gether with the generally sub-solar metallicities measured 
in GRB hosts casts serious doubt on the direct relation be- 
tween GRBs and star formation, which is dominated by 
stella r formation in massive (1 < z < 3, see Caputi et al. 
20051) or in interme diate massive (0.4 < z < 1, see 



Hammer et al.ll2005|) galaxies. However, as discussed be- 



low, our finding of a close spatial association and possibly 
a dynamic connection between GRBs and a "nearby" su- 
per star cluster or a massive star-forming region, suggests 
a more indirect connection between GRB and star forma- 
tion. 



4.2. Inferences on GRB progenitors 

What mechanisms can explain these powerful events ? 
Recall tha t the coUapsar mode l favours a WR origin for 
the GRB ijWooslevet al.lll999l: iMeszaroslEool . This is 
supported by the specific nature of the associated super- 
novae, whose type (SNIb or c) indicates that they have lost 
their hydrogen and/or helium envel ope, as is expected if 
they originate from WR stars (e.s. ' Hirschi et alJ bOOfil) . 
This is consistent with the strong Nitrogen excess we ob- 



serve in the SN1988bw remnant region, which may be pro- 
duced by the progenitor wind or the SN of an initially fast 
rotating star (see Mevnet & Maedcr 2002). 

Are WR stars found in all GRB host galaxies ? We 
have investigated archival data from ESO telescopes, to 
acquire good spectra of the most nearby GRB host galax- 
ies, i.e. those for which 8 meter class telescopes are 
able to detect the faint signature of this stellar popu- 
lati on. Considering t he four similar objects pointed out 
by ISoderberg et alJ l|2004hl see their Fig. 2), this in- 
cludes our two targets plus GRB031203 (z = 0.1055) and 
GRB030329 {z = 0.169). Altogether these objects repre- 
sent the low energetic end of GRBs, and for all of them, 
supernovae have been observed few weeks after the GRB. 
Nevertheless, GRB020903, GRB030329, and GRB031203 
are almost 100 times brighter than GRB980425 and are 
intermediate between this event and cosmological GRBs. 
Unfortunately available VLT data on the GRB030329 host 
are of insufficient quality ^ . On the other hand, the spec- 
trum of the GRB031203 host has a quality similar to our 
spectrum of GRB020903. Figure 5 shows its spectrum and 
we can see that a blue bump around the Heii4686 line and 



^ We extracted deep observations made at VLT in excellent 
seeing conditions (0.5 arcsec); unfortunately this very compact 
object has been observed with a wide slit (1.3 arcsec) and the 
region near the bump is highly contaminated by residuals of 
the strong Ol5577 sky line. 
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also a bump around the Civ5808 line are unambiguously 
detected, indicating the presence of WR features, respec- 
tively of the WNL and WC stars. For all the observed 
GRB hosts having a spectrum with enough quality to de- 
tect WR stars, we indeed detect WR stars. 

Our observations of the WR region of the GRB980425 
host show a very compact (5 pc FWHM) super star cluster 
containing thousands of O stars and almost one hundred 
WR stars. The bright region in GRB020903 separated by 
a projected distance of ~ 460 pc from the SN (see Fig. 4) 
is probably responsible for at least a part of the observed 
WR emission. It also includes a compact component, ap- 
parently not resolved by the ACS camera and therefore 
smaller than 100 pc at z = 0.25. Unfortunately there is 
no HST imagery of the GRB031203 host, and the precise 
location of the GRB/SN is still unclear from ground-based 
data (GRB from 0.2 to 1 arcsec off the host galaxy center 
according to Gam-Yan et al. 2004 or to P rochaska et al. 
2004) . Images of the GRB030329 host (see lFruchter et all 
l2003j) reveal a very faint host galaxy {My = —16.5) and 
the SN is found at its edge, off by 750 pc from the galaxy 
center (see Fig. 4). The offsets revealed by Fig. 1 and 4 are 
very robust because they are based on HST images taken 
at different epochs revealing the SN, and the accuracy is 
better than half an ACS pixel ( 0.05 arcsec). They d o not 
differ much from that found bv lBloom et alJ <|2n02|) on a 
larger sample of GRBs (median value of 1.3 kpc for the 
offset), although their result is less accurate because it is 
based on a comparison between HST and ground based 
images. Thus it appears that GRBs often (or always ?) 
occur in the outskirts of regions populated by hot and 
massive stars including WR stars. 

Based on this new information on the spatial lo- 
cation and the stellar populations we now discuss the 
possible scenarios for GRB progenitors. Assuming that 
GRB980425 is a prototype of other GRBs discussed in 
this paper, we are left with two hypotheses: either the 
GRB progenitor is born in a low stellar density region (in 
situ hypothesis) or it has been expelled from high density 
stellar regions (runaway hypothesis). 



4.3. GRB as runaway, fast rotating massive stars 
expelled from superstellar clusters 7 

If GRB hosts include clusters with WR stars in large num- 
bers, why does the SN (and then the GRB) occur several 
hundreds of parsec away from the WR region ? Assuming 
this is not a coincidence we suggest the following "runaway 
GRB" scenario, which may be essential to achieve the nec- 
essary conditions for the formation of a GRB. Within the 
high stellar density of massive super star clusters, some 
stars (or double stars) can be ejected dynamically after 
one or more elastic collision or from supernova kicks in bi- 
nary systems. During such a collision, the star may have 
acquired a very large angular momentum, enough to lose 
its hydrogen or even helium envelope providing, typically 



~ 3 Myr later, a progenitor of a SNIbc such as those ob- 
served after the GRB event. 

In the case of GRB980425, for the progenitor to travel 
800 pc from the WR region to the SN1998bw region in 
say 3-6 Myr, a velocity of ^ 260-130 km s^^ is neces- 
sary. Such velocities are somewhat larger than (but not 
unseen) the typical velocity of Galactic O to early B-type 
runaway stars, which are thought to be runaway stars from 
stellar clusters (see e.g. lBlaauwlll993l: iTenies et al.ll200ll: 
Ide Wit et alJl2005^ . This leaves this possibility for suffi- 
ciently young, i.e. massive GRB progenitors. 

We do not suggest that all massive runaway stars give 
rise to a GRB, as this would correspond to ~ 10-15% of all 
O stars, based on the knowledg e of Galactic runaway stars 
l|Giesll987trde Wit et al.l2005>l . However, in the case of the 
very massive and compact super star clusters observed in 
the GRB host galaxies, the dynamical conditions, ejec- 
tion probabilities, and the resulting properties of runaway 
stars may be quite different from the less dense and less 
populated assoc iations and clusters typically found in our 
Galaxy (see e.g.lLe onard & Du ncanlll99nl : iLeonard Ill99li 
IPortegies Zwart et al...l999.') . Furthermore, only the run- 
away stars with some peculiar ejection history and high 
angular momentum may produce a GRB. An interest- 
ing case possibly resembling our suggested runaway stars 
is the high velocity star HIP 60350 thought to be dy- 
namically eje cted from the compact region NGC 3603 
IjTenies et al.M2001 '). However, this object is of spectral 
type B3-4V, i.e. a priori of too low mass t o produce a 
collap sar via "normal" single star evolution. Hob bs et alJ 
1I2OO5I) also find large velocities for the proper motions of 
233 pulsars with mean ID speeds of 152 km s^^ {a ~ 265 
km s~^). 

Runaway stars can be produced either by supernova 
explosions in massive close binary system or vi a stron; 
dynamical interactions in young star clusters I Blaau' 



w| 



■1993.') . It is beyond the scop e of this paper to differ e ntiate 
between these alternatives. IPortegies Zwart et al.l l|1999|) 
have simulated the "ecology" of a similar system, the cen- 
tral R136 cluster of 30 Doradus in the Large Magellanic 
Cloud. They fomid that physical collisions between stars 
are quite frequent, and closely linked with the evolution 
of the star cluster. Portegies Zwart et al. (1999, see also 
Leonard 1995) argue that dynamically ejected runaway 
stars can be massive and could have acquired a large an- 
gular momentum which might be consistent with expec- 
tations for a SNIbc progenitor. iBallv fc Zinnecken Ij2005^ 
suggest that the merger of two massive stars within dense 
clusters could be a pathway leading to hypernovae/GRBs. 
Combining our empirical findings with their models sug- 
gests that GRBs originate from (very) rapidly rotating 
stars - resulting from previous stellar collisions in the 
cluster core - which are ejected during a subsequent dy- 
namical interaction. Such dynamical interactions might 
also help to alleviate the difficulties to form GRBs in ro- 
tating stellar models including the breaking by magnetic 
fields found in particular a t solar o r somewh at subsolar 
metallicities IjPetrovic et al.ii2005: Wooslev fc Heger 20061 
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lYoon k, Lange]il2005|) . Most likely the runaway GRB pro- 
genitor ejected from the cluster is a single star or a tight 
binary. 



This scenario has several advantages. It explains (by 
construction) the observed spatial shift between the GRB 
position and a nearby massive star-forming region, as 
observed for GRB980425 and also tentatively indicated 
for GRB030329 and GRB020903. It accounts for the fact 
that all our studied GRB host galaxies are indeed WR 
galaxies. Furthermore, it allows us to reconcile single 
star collapsar models with the lack of massive stars ob- 
served in the immediate vicinity of GRB980425. Indeed, 
a rather low branching ratio of GRB/SNIb c (typically 
R = j y(GRBSNe)/JV(SNIbc) ~ (2-4) x IQ-^. lvan Putl^ 
I2OO4I: IPodsiadlowski et al.ll2004|) implies that statistically 
a GRB should be accompanied by a large number of mas- 
sive stars, which are not found in this region. In other 
words, if GRBs result from the tail of a distribution of 
properties of massive single stars, the population corre- 
sponding to the remainder of this distribution, i.e. sev- 
eral thousand massive stars, should be present. The age 
and the metallicity of the super star cluster (SSC) is com- 
patible with massi ve single star progenitor models (e.g. 
iHirschi et ^]Em^. 



On the other hand, if the progenitor of SN1998bw was 
ejected from the nearby SSC ("WR region" in Fig. 1) in 
a random direction, why is it observed within a small, 
but relatively inconspic uous region surroun ded by 6 point 
sources (see Fig. 1 and lFvnbo et al.ll2000l) ? Using Fig. 1 
we have tested if this spatial association is purely fortu- 
itous, assuming that the GRB progenitor is a runaway 
star expelled from the WR region. We have drawn a cir- 
cular annulus centered on the WR region with an internal 
radius of 0.55 arcsec (to exclude the whole WR region) 
and an external radius of 5.3 arcsec (to include the SN re- 
gion); it corresponds to an area of 87.3 arcsec^. Assuming 
that the runaway star has been expelled in a random di- 
rection, we calculate the probability of the association of 
GRB980425 with the SN region. Wc find that 9.6 arcsec^ 
among 87.3 arcsec^ correspond to pixels with larger sur- 
face brightness than the SN region, i.e. to denser stellar 
regions. It leads to a marginal probability (11%) of finding 
an expelled runaway star in a region as bright or brighter 
than the SN region. A similar but more complex estimate 
can be attempted by accounting for all the point sources 
found in the corresponding area. It is, however, limited by 
several bright regions which could not be resolved into in- 
dividual sources. The number of point sources in the 87.3 
arcsec^ area ranges from 250 and 300. The probability of 
finding the expelled star surrounded by 6 point sources 
within 1 arcsec^ ranges from 7 to 13%. We conclude that 
the location of the SN is marginally fortuitous and thus 
not unexpected. 



4.4. GRBs from binaries occurring in situ within low 
stellar density regions ? 

If the GRB progenitor was born in a low density region, 
it is rather implausible that it is a single, fast rotating, 
massive star. Alternatively, several binary scenarios have 
been proposed for the form ation of long duration GRBs 
fe.g. lFrver fc Wooslevlll999t) . Ages of ~ 6-8 Myr (see sec- 
tion 3.1.2) are a priori compa tible w ith binary collapsar 
or He star merger models f see iFrver & Wooslcv 1999). 

A (massive) binary scenario does not require a very 
massive, rich cluster or region, as the low GRB/SNIbc 
branching ratio might e.g. be explain ed by processes re - 
lated to the nucleation of black holes Ijvan PutterJl2004f) . 
However, there is no reason why in this case GRBs should 
be found preferentially in low stellar density regions (e.g. 
less than 10 O stars in the SN region of GRB980425). 
When considering larger samples, one should then find 
GRBs spatially distributed proportionally to the number 
of massive stars, i.e. to current star formation. 



4.5. Further observational tests 

Our observations and the two scenarios discussed above 
imply that the "classical" single star collapsar model has 
to be abandoned for a binary scenario or for our newly 
suggested "runaway ejection scenario" (or a combination 
of both) . What additional and possibly decisive tests dis- 
tinguishing these scenarios can be envisaged? 

The most decisive test of the runaway ejection scenario 
would be the localisation of more isolated GRBs together 
with a "parent" super star cluster from which they were 
ejected. The lack of a fairly massive and young cluster in 
reasonable proximity of the GRB would exclude our sce- 
nario. We are also puzzled by the fact that, at least in 
two cases (GRB 980425 and GRB030329, see Figure 4), 
the GRB/SNs seem to be located in the prolongation of 
the elongated WR regions. Is this consistent with a run- 
away scenario ? Such investigations require the detection 
of many nearby GRBs, preferentially at z< 0.05, which is 
challenging. 

Other potentially interesting constraints may 
come from the study of absorpt ion lines originating 
from the burst env i ronment (cf. 

'Mirabal et al." '2003'; Starl ing et all , 

2006: Prochaska et al. 20061) . For instance, narrow Fe li* 
and Si il* fine structure absorption lines have recently 
been observed in GRB051111, indicative of a high density 
medium in close proximity to a SSC or to the GRB 
i Berger et aLll2006t |p7ochaska et alll2006^ . ISerger et alJ 
2006j argue that the source of radiative pumping 
required for the excitation of the fine structure levels 
must be located close to the absorbing medium on the 
GRB line of sight. If this exciting source is a SSC, as 
they propose, the GRB is supposed to be located within 
the cluster - in contrast with our observations -, or there 
must be a chance alignment of the GRB and the SSC 
with its absorbing medium, which seems quite unlikely. 



Schaefer et all l2003l:_ 
'20051 iBerger et alJ 
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Alternatively, IProchaska et al.1 l)2006|) argue that the 
absorbing medium could be circumstellar gas around the 
GRB progenitor star. Maintaining such a high density 
circumstellar envelope around a single star ejected with 
high velocity seems quite improbable. If such features 
turned out to be ubiquitous, our suggested scenario would 
probably require the ejection of a tight binary rather 
than that of a single star. 

Clearly more detailed observations including accurate 
spatial information are needed to progress on these issues. 

5. Conclusion 

We have reported here the detailed properties of the en- 
vironment of the four closest GRBs known to date, all 
of them being associated with SNe. For three of them 
(GRB980425, GRB020903, and GRB031203) we have 
spectra of sufficient S/N allowing us to examine the pres- 
ence of WR stars, and all of them show the character- 
istic WR bump. HST images are also available for the 
GRB980425, GRB020903 and GRB030329 hosts. 

WR galaxies arc thought to be related to very re- 
cent starbursts (< 10 Myr). It is a very heteroge- 
neous class of galaxies, o ccurring in var ious en vironments 
i Schaerer fc Vaccall998l) . Following the lContil |l991l) clas- 
sification (more than 100 WR stars per W R galaxy) the 
GRB hosts studied here are WR galaxies. ISmithI l)l99l[l 
estimated that 10 % of HII galaxies are indeed WR galax- 
ies. This somewhat supports the link between GRB and 
WR stars, i.e. GRBs with very recent formation of massive 
stars. 

However, GRBs occurred at several hundreds of parsec 
away from the massive star forming regions hosting large 
numbers of O and WR stars. We notice that this applies 
to all the GRB hosts for which we have enough spatial res- 
olution (see Fig. 4). Assuming this is purely a coincidence 
is not very plausible and would leave open the question 
of the formation of GRB/SNIbc in a medium containing 
only a small number of massive stars (e.g. less than 10 O 
stars within a radius of 100 pc around the SN1998bw re- 
gion) . We believe that such a coincidence is not fortuitous 
and propose therefore a new scenario, where GRBs are as- 
sociated with runaway stars ejected from regions showing 
strong and recent star formation. 

Our observations support the coUapsar model by re- 
vealing for the first time the link between WR stars and 
GRBs associated with SNIbc, although in a less direct 
way than could have been expected. In our proposed sce- 
nario the frequency/rareness of the GRB events could be 
related to dynamical processes occurring during the clus- 
ter ejection. It is also very likely that additional parame- 
ters, such as metallicity and others suggested earlier, play 
a role in establishing the necessary conditions for a GRB 
event. The location of GRBs in sub-lumino us galaxies and 
the absence of GRBs in massive galaxies tLe Floc'h et al.l 
l2003ll2006HFruchter et"allll999|) . as well as the measure- 
ment of sub-solar metallicities in GRB hosts and nearby 
regions (e.g. iProchaska et al. 2004 and this paper), show 



that metallicity also plays a role. This question will be 
addressed in more detail in a subsequent paper. 

The discovery of a possible link between runaway WR 
massive stars and GRBs requires further confirmation. 
Massive binary scenarios must also be explored and tested. 
It has to be understood if such mechanisms occur prefer- 
entially in regions with a low density of stars or not. We 
hope that further nearby GRBs will be soon discovered by 
Swift and other satellites. Follow-up should include very 
high resolution images from the HST at several different 
epochs, and deep spectroscopy of the region hosting the 
GRB. X-shooter at VLT, equipped with an integral field 
unit will be ideal for the latter purpose. 

Our results illustrate the need for a better understand- 
ing of the physics of massive stars and the importance of 
dynamical processes in cluster environments. It is still un- 
clear whether the above results support the use of GRBs 
to test cosmological models. 
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